Transition region properties of a trapped quasi-two-dimensional degenerate Bose gas 
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The classical field simulation technique is used to study a trapped quasi-two dimensional Bose gas. We 
calculate the central curvature of the system density and fluctuations of the condensate mode in the degenerate 
regime. These results provide new understanding of the system behavior in the region of the superfiuid transition. 
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I. INTRODUCTION 

Evidence for the Berezinskii-Kosterlitz-Thouless (BKT) 
superfiuid transition in a dilute quasi-two dimensional (quasi- 
2D) Bose gas was first reported by the ENS group [1 1. In that 
work the algebraic decay in first order correlations, inferred 
from an interference procedure, was used to identify the tran- 
sition point. Subsequent work suggested that the transition 
could also be identified by bimodality in the system density 
after expansion [2j. In contrast, recent experiments at NIST 
[3 1 found that the onset of bimodality and superfluidity are 
distinct, with bimodality occurring at a higher temperature. 

Meanfield methods are inapplicable to the regime of inter- 
est and theoretical understanding of this system is just begin- 
ning to emerge with the development of classical field (c-field) 
[4 1 and quantum Monte Carlo (QMC) 1 5 1 methods for simulat- 
ing the trapped quasi-2D Bose gas. Previous results obtained 
with the c-field technique have found the onset of bimodal- 
ity and BKT superfluidity to be distinct, in qualitative agree- 
ment with the NIST observations |6|. We also note work on a 
semiclassical field method |7|, which was applied to studying 
thermal vortices in a homogeneous 2D Bose gas. 

Many important questions remain to be answered for a 
complete understanding of this system, particularly in the 
transition region. Here we address two important issues: First, 
we examine the use of the central density curvature (n) as a 
means to identify the transition point. This quantity was orig- 
inally proposed by Holzmann and Rrauth who showed, 
using QMC calculations, that n was peaked in the vicinity 
of the superfiuid transition. However, their study contained 
only a small number of results across the transition region, 
and was unable to resolve any distinction between bimodality 
and the onset of BKT superfluidity. The efficiency of the c- 
field method allows us to characterize the behavior of k with 
fine resolution across the transition region. In addition to pro- 
viding us with a better understanding of this observable in the 
region of current interest, it also provides a useful qualitative 
comparison of the c-field and QMC [5 1 techniques. Second, 
we present results characterizing the properties of the system 
in the temperature range where the gas is bimodal but has 
not yet attained BKT superfluidity. In previous work [6| we 
showed that the system has a small but finite condensate frac- 
tion in this temperature range. However, in this regime the 
condensate appeared to fluctuate strongly and easily be pene- 
trated by single vortices (also see 0). We apply a newly de- 



veloped technique E) to sample the condensate number dis- 
tribution and reveal its fluctuations. These new results shed 
important light on the physics of the trapped quasi-two di- 
mensional Bose gas, and suggest some new directions for ex- 
perimental investigation. 



H. FORMALISM 

An ultra-cold Bose gas is described by the Hamiltonian 

H = J d 3 x^( x ) j# sp + ?!^!L$t( x )$( x ) j $( x ) ( 

(1) 

where m is the atomic mass, a is the s-wave scattering length 
and the single-particle Hamiltonian is given by 

u P 2 ,l/22i 2 2, 2 2\ /o\ 

H s P = ^ + -m(uj x x +u y y +uj z z). (2) 

In the regime where the trap frequencies satisfy lu z 3> ui x ,u> y , 
motion in the z-direction may be frozen out, resulting in a 
quasi-2D regime. Eq. ([T]) assumes three-dimensional scatter- 
ing between the atoms which requires that the z -confinement 
length, a z = ^h/muj z , satisfies a z 3> a (9). These require- 
ments are well-satisfied by experiments lITl l2l l3l [TOl . 

Our simulation method is based on a c-field representation 
of the highly occupied low energy modes (i.e. those below an 
appropriately chosen energy cutoff, see IfTTI ) and a Hartree- 
Fock treatment of the remaining sparsely occupied (high en- 
ergy) modes of the system. An advantage of this method is 
that the critical modes are contained within the c-field descrip- 
tion and are treated non-perturbatively. For details of this c- 
field theory refer to |[T2"1 [T3~l H4l and for the specific applica- 
tion to the quasi-2D trapped Bose gas see Bl l6l [T5l [161 . 



III. RESULTS 

We simulate 87 Rb atoms in a cylindrically symmetric trap 
(i.e. Lo x y = ui), presenting results for two sets of parame- 
ters in which the dimensionless 2D interaction strength, g = 
\/&Ka/a z 13, takes the values 0.075 and 0.107. For com- 
parison, the ENS experiment (with 87 Rb) was in the regime 
where g w 0.13 HJ and for the NIST experiment (with 23 Na) 
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Figure 1: Fraction of atoms in the quasicondensate (dots), con- 
densate (squares) and superfluid (triangles) components as a func- 
tion of reduced temperature. Figure from Ref. [6]. Parameters: 

(Lo XyV ,uj z ) = 2tt x (9.4,1880) Hz, N = 13.5 x 10 3 atoms at 



Tbkt = 34.7 nK, with i 



0.107. 



g w 0.02 0. In this paper the densities given are areal (i.e. in- 
tegrated over the z-direction) and are dependent on the 2D 
position vector r = (x, y). 



A. Degenerate components 

In previous work 0] we identified the various components 
contributing to the phase diagram. Here we briefly review the 
definitions of these components and summarize those results. 

a. Quasicondensate is the component of the system 
with suppressed density fluctuations, defined by ifTTl 



n qc (r) = V2(n(r)) 2 -(n(r)2), 



(3) 



where n(r) is the density operator for the system. Note, a 
normal system with Gaussian fluctuations has n qc = 0. 

b. Condensate/Coherence is identified by a macro- 
scopic eigenvalue in the one-body density matrix G(r,r') = 
(^(r)^! - ')), ||5][T3]. The condensate is expected to vanish 
for the uniform 2D gas in the thermodynamic limit, however 
our results show that it plays a role in the finite trapped sam- 
ples realized in experiments (also see lB1 [l8l ). 

c. Superfluid model of Holzmann and Krauth In Ref. 
a model for the superfluid component was proposed, based on 
a local density application of the uniform results IfTTl [T9l to 
the trapped system. This model predicts a superfluid compo- 
nent wherever the local density exceeds the critical value, 



= A- 2 log 



(4) 



with C = 380 ± 3 and A = h/y/2Trmk B T flj). The tempera- 
ture at which the peak phase space density satisfies the critical 
condition Q is denoted Tbkt lf51 l6l l20l |2T1 1 . 




Figure 2: Total areal density profiles at several temperatures: n(x) 
(solid curve), fitting region (vertical lines indicate x = ±5/i m) and 
quadratic fits (dashed). Parameters as in Fig.[T| 



The results of Ref. [6 | for these components are summa- 
rized in Fig. [T] plotted against reduced temperature, T' — 
T/Tq, where T is the ideal gas condensation temperature 2 . 

These results clearly reveal that the quasicondensate, con- 
densate and superfluid components are distinct in the quasi- 
2D system. Furthermore, we observe that, in order of decreas- 
ing temperature: (i) The first manifestation of degeneracy is 
that density fluctuations are suppressed, i.e. a quasiconden- 
sate forms (this begins to occur at T' > 1 in Fig. [TJ. (ii) At 
T c a condensate appears and the momentum distribution of 
the system becomes bimodal. This crossover is a finite size 
effect, but occurs well-separated from the BKT transition for 
the typical experimental regime. We identify T c as where the 
largest eigenvalue of G(r,r') exceeds the second largest by 
a factor of 1.5, at which point momentum space bimodality 
is clearly apparent |6). (iii) At Tbkt the peak phase space 
density is sufficiently high for BKT superfluidity to emerge in 
the system. More details on these results and the simulation 
method are given in Ref. J6) . 



B. Central curvature of density distribution 

Holzmann and Krauth |5| proposed using the central den- 
sity curvature 



dn(x) 
' d(x 2 ) 



(5) 



as a signature of the superfluid transition, showing with QMC 
calculations that n was peaked near the superfluid transition. 
However they only calculated a small number of results across 
the transition region and it was not possible to clarify if this 
peak occurred at T c or Tbkt- 

In Fig. [2] we show two examples of our c-field calculated 
average density profiles, and the quadratic curves fitted to the 
central region which we use to extract the curvature. In gen- 
eral we observe that the quality of the fits become progres- 
sively poorer with increasing temperature. We find that the 



1 We use the lowest axial mode areal density to identify Tbkt, whereas 
1201 use the total areal density, also see 1211 . 



2 An improvement in our meanfield treatment of the incoherent region atoms 
has led to a small change in the fractions compared to those reported in (5). 
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Figure 3: Central density curvature. Parameters: (a) (tJz,,,, tj z ) = 
2tt x (9.4, 1880) Hz with g = 0.107, (black dots) system with N = 



21.2 x 10 3 at Tbkt = 43.1 nK. (b) (a 



= 2tt x (9.4, 940) 



Hz with g = 0.075, (black dots) system with N = 30.2 x 10 d at 
Tbat = 46.8 nK, (gray dots) system with N = 67.3 x 10 3 at 
Tbkt = 64.9 nK. Vertical solid and dashed lines indicate Tbkt 
and T c , respectively. The Horizontal dashed line shows ktf- 



curvature evaluated this way depends on the size of the fitting 
region, and we use \x\ < 5/j,m. 

Figure[3]presents results for the central curvature as a func- 
tion of temperature. We find, in qualitative agreement with 
the QMC results in Ref. [5], that the curvature is peaked in 
the transition region. For low temperatures (T < Tbkt) 
the curvature approaches a constant at a value slightly greater 
than that of the Thomas-Fermi profile (ktf = m 2 uj 2 /2h 2 g). 
This is consistent with the system density fluctuations being 
strongly suppressed at trap center. Above the transition region 
(T ~ Tbkt) we observe the curvature decrease with increas- 
ing temperature. 

The range of results in Fig. [3] span a considerable range 
of system parameters, and show that the curvature peak near 
Tbkt is robust. Our results also show that as the system size 
increases Tbkt and T c gradually get closer (e.g. compare the 
black and gray results in Fig. [3]). However, our results do not 
allow us to conclude whether this is because of a reduction in 
the finite size effect (due to larger atom number) or because 
the system is crossing over to being three-dimensional (due to 
the increase in transition temperature relative to hutz/ks)- 
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Figure 4: Condensate number distribution: Histograms of the con- 
densate number extracted from c-field evolutions at a range of tem- 
peratures. The " x4" indicates that the frequencies have been divided 
by 4 for subplot (h) for convenience. Parameters as in Fig.[T] 



ical exponents and emergence of coherence for a three- 
dimensional trapped Bose gas. Here we apply these tech- 
niques to the quasi-2D Bose gas. 

The procedure for extracting the condensate number statis- 
tics is discussed in detail in [8|. Briefly, we obtain samples 
of the condensate amplitude (a c (t)) at time t, by taking the 
inner product of the classical field at this time with the (time- 
averaged) condensate mode 3 . We can then construct the con- 
densate number distribution, with number samples given by 
N c (t) — |a c (i)| 2 . For each simulation we evaluate 5 x 10 4 
samples of N c taken over 1000 radial trap periods of c-field 
evolution. Figure[4]shows the resulting histograms of the con- 
densate number for eight different temperature regimes. 

The results at the lowest [Fig. |4ja)] and highest [Fig. |4jh)] 
temperatures show number distributions consistent with co- 
herent and incoherent number statistics, respectively. What 
is most intriguing is the qualitative change in behavior we 
observe in the temperature range from Tbkt to T c . For 
T ~ Tbkt [Fig. |4j e )l we find tnat tne distribution is max- 
imum at a macroscopic number (N c max w 700), with negli- 
gibly small probability of finding iV c = 0. For temperatures 
between Tbkt and T c the condensate number distribution has 
an interesting shape: the most likely condensate number is a 
macroscopic, non-zero value, but the probability of finding 
N c = is appreciable (e.g. see Figs.|4|f) and|4^g)). In this 
regime condensate number fluctuations are large relative to 
the mean value. Figure|5]demonstrates a N c (t) trajectory from 



Condensate number fluctuations 



Recently the c-field technique has been extended to probe 
mode statistics in detail, and applied to exploring the crit- 



' Taken as the eigenvector of the one-body density matrix corresponding to 
the largest eigenvalue, which above T c approximately corresponds to the 
harmonic oscillator ground mode. 
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Figure 5: (upper) Condensate number variation during c-field evo- 
lution for system at T' — 0.91, i.e. same data as for Fig. |4jf). 
(lower) C-field density of the A,B microstates at the times indicated 
in the upper plot. Inner (outer) circle indicates the condensate (quasi- 
condensate) 1/e radius. Vortices (+) and antivortices (-) indicated 
inside the central condensate region only. 

a c-field simulation for this regime. This trajectory reveals a 
typical event that occurs at this temperature: a sustained pe- 
riod (50.2 s < t < 50.7 s) over which iV c is suppressed to a 
value significantly below the long time average. Examining 
the microstates we find that this suppression occurs when an 
isolated vortex enters the central region of the system (see mi- 
crostate B in Fig. [5]). These free vortices are long lived and 
for Tbkt < T < T c they penetrate the central region quite 
frequently. Below Tbkt such free vortices are strongly sup- 
pressed (near trap center) leading to a negligible probability 
of observing N c « (of course this does not exclude paired 



vortices). The role of vortices was speculated to be the origin 
of fluctuations observed in the NIST experiment [3] for tem- 
perature above Tbkt, although limited resolution prevented 
the direct observation of these vortices (also see| 10|). 



IV. CONCLUSIONS 

We have analyzed two important characteristics of the low- 
temperature trapped quasi-2D Bose gas. Our results for the 
central density curvature of the system show that this quan- 
tity peaks near the transition region, in qualitative agreement 
with previous quantum Monte Carlo work. The feature of our 
work is that our results resolve the transition region in detail 
and clearly reveal that the curvature peak occurs in the vicinity 
of Tbkt (and not T c ), and is thus distinct from the momen- 
tum space bimodality arising at T c . So far most experiments 
have examined the system density profile after expansion for 
intermediate times, where the imaged density is a convolu- 
tion of the in situ position and momentum space distributions 
(e.g. see 11221 ). This suggests that the use of longer time-of- 
flights (or Bose-gas focusing [23 1) to measure the momentum 
distribution, and in situ imaging to directly observe the posi- 
tion density, will more clearly reveal the distinction between 
T c and Tbkt in the quasi-two dimensional gas. 

Our results for the condensate fluctuations help clarify the 
changes that occur when coherence and superfluidity emerge 
in the system. In particular, it reveals that the condensate 
above Tbkt is frequently penetrated by free vortices causing 
strong number fluctuations. These predictions broadly support 
the interpretation presented in Ref. [3|. To date no direct ex- 
perimental measurements of the condensate number statistics 
have been made. However, with the emergence of new tech- 
niques for accurately measuring critical properties [24], and 
proposals for the use of related observables (e.g. fluctuations 
in the zero momentum density J8)) to infer this information, 
we expect such measurements will be feasible in the future. 
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